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Abstract 



In this paper we present the results of a search for nuclearites in the penetrating 
cosmic radiation using the scintillator and track-etch subdetectors of the MACRO 
apparatus. The analyses cover the f3 = v/c range at the detector depth (3700 
hg/cm2) 10^5 < /5 < 1; for /? = 2 X lO"^ ^he flux limit is 2.7 x 10"^*^ cm-^s-^sr"! 
for an isotropic flux of nuclearites, and twice this value for a flux of downgoing 
nuclearites. 



1 Introduction 



In 1984 Witten formulated the hypothesis that "strange quark matter" 
(SQM) composed of comparable amounts of u, d and s quarks might be the ground 
state of hadronic matter. "Bags" of SQM (also known as "strangelets" ) would be 
heavier than a single A'^ baryon, but neutral enough not to be limited in size by 
Coulomb repulsion. 

If particles of SQM were produced in a first-order phase transition in the early 
universe, they would be candidates for Dark Matter (DM), and might be found 
in the cosmic radiation reaching the Earth. Particles of SQM in cosmic radiation 
are commonly known as "nuclearites" 0. This paper describes an experimental 
search for nuclearites with the MACRO detector at Gran Sasso. 

The main energy loss mechanism for nuclearites passing through matter is 
elastic or quasi-elastic atomic collisions The energy loss rate is 

where a is the nuclearite cross section, v its velocity and p the mass density of the 
traversed medium. 

For nuclearites with masses M > 8.4x10^^ GeV/c^ (~ 1.5 ng) the cross section 
may be approximated as 

/ 3M \2/3 

cr ~ vr X (2) 



V47rp 



'N 



where pN (the density of strange quark matter) is estimated to be pN — 3.5 x 10^^ 
g/cm^, somewhat larger than that of atomic nuclei 0]. For nuclearites of masses 
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M < 8.4 X 10^^ GeV/c^ the collisions are governed by their electronic clouds, 
yielding o" ~ vr x 10~^^ cm^. 

An experimental search for nuclearites has an acceptance that depends on 
nuclearite mass, since upgoing nuclearites that traverse the Earth before detection 
only occur for sufficiently large nuclearite masses ( ~ 6 X 10^2 GeV/c2 at typical 
galactic velocities Q). 

An upper limit on the nuclearite flux may be estimated assuming that 
^max. = Pdmv /{2'kM), where pdm — 10^^^ g/crc? represents the local 
DM density, and M and v are the mass and the velocity of nuclearites, 
respectively. 

Different indirect methods to search for nuclearites have been suggested 
[|]. Some exotic cosmic ray events were interpreted as due to incident 
nuclearites, for example the "Centauro" events and particles with anomalous 
charge/mass ratio [4-10]. The interpretation of those possible signals is not 
yet clear. Searches for strangelets are being performed at the Brookhaven 



AGS fill and at the CERN-SPS fl2 



Relevant direct flux upper limits for nuclearites came from two large 
area experiments using CR39 nuclear track detectors; one experiment was 
performed at mountain altitude |jl3|, the other at a depth of 10"^ g/cm^ at 
the Ohya stone quarries |I4 |. 

The lowest flux limits have been obtained by examining ancient mica 
samples. It should be kept in mind however that this technique has inherent 
uncertainties [|15|, |16]. 

MACRO (Monopole, Astrophysics and Cosmic Ray Observatory) is an 
underground detector located at the Gran Sasso Laboratory in Italy, at an 
average depth of 3700 hg/cm^ and at a minimum depth of 3150 hg/cm^. 
MACRO uses three different types of detectors: liquid scintillators, limited 
streamer tubes and nuclear track detectors (CR39 and Lexan) arranged in 
a modular structure of six "supermodules" (SM's). The overall dimensions 
of the apparatus are 76.5 x 12 x 9.3 m^ [0. The response of the three 
types of detectors to slow and fast particles has been experimentally studied 
[18-20]. One of the primary aims of MACRO is the search for superheavy 
GUT magnetic monopoles |2T|, E^. Some of the search methods used for 



this purpose, namely those based on the liquid scintillators and nuclear track 
detectors, may also be applied to search for nuclearites. 

In computing the acceptances for nuclearites of the scintillator and of the 
nuclear track subdetectors one has to take into consideration the nuclearite 
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absorption in the Earth. In Fig. 1 we present the fraction of sohd angle from 
which nuclearites might reach MACRO as a function of their mass, assuming 
a velocity at the ground level of /3 = 2 x 10~^. For masses smaller than 5 x 10^^ 
GeV/c^ the nuclearites cannot reach the detector; for 5 x 10^^ < M < 10^^ 
GeV/c^ only downward going nuclearites can reach it; for M > 10^^ GeV/c^ 
nuclearites from all directions can reach MACRO. The detector acceptance 
for an isotropic flux (at ground level) of nuclearites has to be scaled according 
to such curves. 

Some earlier results obtained using the scintillator subsystem of the lower 
part of the first SM have been published |]23|. More recent results obtained 



with the scintillator and CR39 subdetectors were reported in Ref. p4 . 

As recently suggested by Kusenko et al. the MACRO search 

for nuclearites could also apply to Supersymmetric Electrically Charged 
Solitons (Q-balls, SECS). Q-balls are supersymmetric coherent states of 
squarks, sleptons and Higgs fields, predicted by minimal supersymmetric 
generalizations of the Standard Model [^; they could be copiously produced 
in the early universe. Relic Q-balls are also candidates for cold DM . 



In the following sections we present results on the search for nuclearites 
in the liquid scintillator and in the nuclear track subdetectors of MACRO. 
The density of the gas in the MACRO streamer tubes [|l9l is too low to allow 
the detection of nuclearites for the hypothesized energy loss mechanism. 

2 Searches using the hquid scintillator 
subdetector 

De Riijula and Glashow have calculated the light yield of nuclearites 
traversing transparent materials on the basis of the black-body radiation 
emitted along the heated track 0]. The light yield per unit track length is 
(in natural units, ft = c = 1) 

% - ^^lZ,i:<nlnrW, (3) 

where m is the mass of a molecule of the traversed material, n is the number 
of submolecular species in a molecule, lOmax is the maximum frequency for 
which the material is transparent, a is the nuclearite cross section and v its 
velocity. 
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In Ref. this formula was applied to our liquid scintillator. It 
was shown that the scintillator subdetector is sensitive even to very small 
nuclearite masses and to low velocities (/? ~ 5x 10^^). The light yield is above 
the 90% trigger efficiency threshold of the MACRO scintillator slow-particle 
trigger system for most nuclearite masses {dL/dX > 10~^ MeV cm~^). 

The scintillators are therefore sensitive not only to galactic {(3 ~ 10"'^) 
or extragalactic nuclearites (higher velocities), but also to those possibly 
trapped in our solar system (/? ~ 10""^). 

The nuclearite detection efficiency in the scintillator subdetector is 
assumed to be similar (or larger) to that for magnetic monopoles; the 
selection criteria used to search for monopole events are also applicable 
for nuclearites. No saturation effects of the detectors, electronics, or 
reconstruction procedure are expected to reduce the detection efficiency of 
the liquid scintillator subdetector. 

Different monopole triggers and analysis procedures have been used in 
the search for cosmic ray strangelets in different velocity domains following 
the evolution of the detector. The relevant parameters of each actual search 
for nuclearites using the scintillator subdetector are presented in Table 1. As 
no candidate satisfied all the requirements, the resulting flux upper limits 
at 90% confidence level (C.L.) are listed in Table 1 and presented in Fig. 2. 
Some details on the present searches are given below. 

The slow monopole trigger which includes the analog Time Over 
Half Maximum (TOHM) electronics and the digital Leaky Integrator 
(LI) electronics [0, recognizes wide pulses or long trains of single 
photoelectrons generated by slow particles, rejecting large and short pulses 
produced by muons or radioactive decay products. When a trigger occurs, 
the wave forms of both the anode and the dynode (for the 1989-91 run 
period) for each photomultiplier tube are separately recorded by two Wave 
Form Digitizers (WFD). A visual scan is then performed on the selected 
events. This procedure was applied to the searches for nuclearites with 
10-5 < /3 < 3.5 X 10-3 H, ||. 



The Fast Monopole Trigger (FMT) is based on the time of flight between 
two layers of scintillators. A slow coincidence between two layers is vetoed 
by a fast coincidence between them. Additional wave form analysis is 
performed on the selected events. This procedure was applied to the search 
for intermediate velocity nuclearites (2.5 x 10-^ < (3 < 1.5 x 10~'^) p3| . 

The scintillator muon trigger. A fast nuclearite should produce a light 
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yield (dL/dx) at least three orders of magnitude larger than that from a 
typical muon. It was checked that no negative effects arise on the detecting 
system from the larger pulse heights. No event was found having a dL/dx in 
both walls greater than 10 times that of a muon. This technique was applied 
in the early analyses for high velocity nuclearites (1.5 x 10~^ < (3 < 1) p3| . 

The Energy Reconstruction Processor (ERP) is a single-counter energy 
threshold trigger |2^, ^ . The ERP analysis requires triggers in two different 
scintillator planes, separated in the vertical direction by at least 2 m, insuring 
a time of flight long enough for accurate velocity measurements. The energy 
deposition must be at least 600 MeV in each counter. The ERP analysis 
was used to search for fast nuclearites (/3 > 0.1) The raw triggering 



efficiency for nuclearites with /3 > 0.1 is essentially 100%. 

Pulse Height Recorder and Synchronous Encoder (PHRASE) is a system 
designed primarily for the detection of supernova neutrinos The 
event selection requires a coincidence between two scintillator planes, with 
no more than 2 contiguous hits in each plane, with an energy release of at least 
10 MeV in each layer. It was checked that no negative effects arise from larger 
pulse heights. A minimum separation of 2 m is required for hits in the two 
counters, while a software cut {f3 < 0.1) is imposed in order to reject the tail 
of the cosmic ray muon distribution. The particle velocity is reconstructed 
using the scintillator time information. The PHRASE search for nuclearites 
covers a large velocity range: 1.2 x 10^^ < (3 < 10^^ [^. The lower limit 



corresponds to the threshold for the detection of bare monopoles with unit 
Dirac magnetic charge {g = gc); as the hght yield produced by nuclearites 
is larger than that of monopoles, the nuclearite search might be extended to 
lower velocities. For candidates with /5 < 5 x 10~^ we compare the duration 
of the scintillation light pulse (measured by the PHRASE WFD) with the 
one computed using the particle velocity; candidates with 5 x 10^^ < /5 < 0.1 
are cross-checked on the basis of the measured energy loss. All the candidates 
with /5 ~ 5 X 10^^ are examined using both techniques, in order to ensure 
the continuity of the analysis. 



7 



3 Searches using the nuclear track 
subdetector 



The nuclear track detector is located horizontally in the middle of the lower 
MACRO structure, on the vertical east wall and on the lower part of the 
vertical north wall. It is organised in modules ("wagons") of ~ 25 x 25 cm^; 
a "wagon" contains three layers of CR39, three layers of lexan and 1 mm 
thick aluminium absorber. Details of the track-etch subdetector are given in 
Ref. the total area is 1263 m^. At the point that this analysis ended 



we had etched 227 m of CR39 with an average exposure time of 7.6 years. 



In Ref. [0 it was shown that the formation of an etchable track in CR39 is 
related to the Restricted Energy Loss (REL) which is the fraction of the total 
energy loss which remains localized in a cylindrical region with about 10 nm 
diameter around the particle trajectory [^]. There are two contributions to 



REL: the electronic energy loss (Se), which represents the energy transferred 
to the electrons, and the nuclear energy loss (Sn), which represents the energy 
transferred to the nuclei in the material. In Ref. it was shown that Sn is 
as effective as Se in producing etchable tracks in our CR39. This result was 
confirmed in Ref. for different types of CR39. In the case of nuclearites 
the REL is practically equal to S^, thus Eq. 1 may be used for calculating 
REL. 

In Fig. 3 we present the energy loss of nuclearites in CR39; the calculation 
assumes that the energy is transferred to the traversed material by displacing 
the matter in the nuclearite path by elastic or quasi-elastic collisions. Such 
processes would produce the breaking of the CR39 polymeric bonds, leaving 
etchable latent tracks, if the energy loss is above the detector threshold. For 
the MACRO CR39 the "intrinsic" threshold is about 20 MeV g-^cm^ in the 
condition of a chemical etching in 8N NaOH water solution at 80° C; this 
is shown in Fig. 3 as the lower horizontal line. The dotted line in Fig. 3 
represents the REL for g = go bare magnetic monopoles in CR39 [^. 

Several "tracks" were observed, mainly due to recoil protons from neutron 
interactions or due to polymerization inhomogeneities. In the conditions of 
the average exposure time in MACRO the number of background tracks is 
about 0.5 /m^ of CR39. When we required that the observed etch cones 
were present on at least four CR39 surfaces and were consistent with being 
from the same particle track, all of the candidates were ruled out. 



8 



From Fig. 3 it is apparent that our CR39 is sensitive to nuclearites of 
any mass and with (5 > 1.5 x 10~^. Nuclearites with mass larger than 
~ 10^^ GeV/c^ can be detected even for velocities as small as /3 = 10~^. 
As a consequence, the 90% C.L. limit for /5 ~ 1 monopoles established by 
the nuclear track subdetector (6.8 x 10"^^ cm~^s^^sr~^) applies also to an 
isotropic flux of M > 5.6 x 10^^ GeV/c^ nuclearites. This limit is presented 
in Fig. 2 as curve "F" and is included in Table 1. For lower mass nuclearites 
the 90% C.L. flux limit is twice this value (1.4 x 10"^^ cm~^s~^sr~^) because 
of the solid angle effect shown in Fig. 1. 



4 Discussions and conclusions 

No nuclearite candidate was found in any of the reported searches. The 90% 
C.L. flux upper limits for an isotropic flux of nuclearites are presented in 
Fig. 2. 

Because either the scintillator or the CR39 can give us a credible 
nuclearite detection, we sum the independent parts of the individual 
exposures to obtain the global limit denoted as "MACRO" in Fig. 2. This 
procedure ensures the 90% C.L. significance of the global limit. 

All limits presented in Fig. 2 refer to the flux of nuclearites at the level of 
the MACRO detector, i.e., below an average rock thickness of 3700 hg/cm^. 
To compare our limit with the limits published by different experiments and 
with the limit calculated from DM density in our galaxy, we integrated the 
energy loss equation for a path corresponding to the averaged rock thickness 
and for different velocities at the detection level. Thus we obtained a relation 
between the nuclearite velocities at the level of the detector and at the ground 
level. Similar calculations were made for other underground experiments 



1-1. l ". 



Fig. 4 shows the 90% C.L. MACRO upper limit for a flux of downgoing 
nuclearites compared with the limits reported in Refs. [jl3[ ( "Nakamura" ) , 
[|g ("Orito"), the indirect mica limits ^ |T§ and with the DM limit, 
assuming a velocity at ground level of (3 = 2 x 10^^. At j3 = 2 x 10~^ 
the 90% C.L. MACRO limit for an isotropic flux of nuclearites is 2.7 x lO"^*' 
cm~^s~^sr~^. In Fig. 4 we extended the MACRO limit above the DM bound, 
in order to show the transition to an isotropic flux for nuclearite masses larger 
than ^ 6 X 10^^ GeV/c^. 
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OIVI a 




Run 
period 


LlVG 

time 

(days) 


iDUtiupic yu /o 

C.L. flux limit 

(cm-2 s-"'^ sr-^) 


Ref. 


A 


TOHM 

FMT 
/i-trigger 


1 


10^5 _2.5 X 10-3 
2.5 X 10-3-1.5 X 10-2 
1.5 X 10-2- 1 


3/89-4/91 


457.5 


5.5x10-15 

4.0x10-1-* 
9.0x10-15 


ii 


B 


TOHM 


1 


10-5 -3.5 X 10-3 


10/89-11/91 


453 


5.6x10-15 


|24 


1 


C 


TOHM 


1-6 


10-5 - 10-2 


12/92-6/93 


163.5 


4.1x10-15 


H 


D 


ERP 


1-6 


0.1 - 1 


12/92-6/93 


166.5 


4.4x10-15 


|24 


1 


E 


PHRASE 


various 


1.2 X 10-3- 0.1 


10/89-12/98 


2587.5 


3.6x10-16 


|24 




F 


CR39 


1-3 


10-5 - 1 


9/88-3/99 


-t 


6.8x10-16 


H* 



* Analyses updated in this paper. 

f For the CR39 it is more appropriate to quote the average exposure time of the etched 
part (7.57 years). 

Table 1. Summary of the nuclearite searches with the scintillator and with 
the CR39 subdetectors. The different techniques are briefly explained in the 
text. Further details are given in the quoted references. At /3 = 2 x IQ-^ 
the combined MACRO limit is 2.7 x IQ-^^ cm-^s'^sr-^. (The procedure 
to obtain the MACRO limit takes care of the overlapping in beta and time 
ranges of individual analyses). 
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Nuclearite mass (GeV/c^) 

Figure 1: Fraction of the solid angle from which nuclearites with 
P = 2 X 10~^ and different masses might reach the MACRO detector. 
For masses smaller than 5 x 10^^ GeV/c^ the nuclearites cannot reach the 
detector; for 10^^ < M < 10^^ GcV/c^ only downward going nuclearites 
can reach it; for M > 10^^ GeV/c^ nuclearites from all directions can reach 
MACRO. 
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10 ^ 10 10 ^ 10 ^ 10 ^ 1 

(S^y/c 

Figure 2: The 90% CL. upper limits for an isotropic flux of nuclearites 
obtained using the liquid scintillator (curves A - E) and the CR39 nuclear 
track (curve F) subdetectors; the bold line is the present MACRO global 
limit. 
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10 10 10 10 

^=v/c 

Figure 3: Restricted energy losses in CR39 for nuclearites of different 
velocities and masses M (in GeV/c^). The dotted line shows the REL of 
magnetic monopoles with unit Dirac magnetic charge {g = go) in CR39 and 
is included for comparison. The threshold of our CR39 is also presented. 
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Figure 4: 90% C.L. flux upper limits versus mass for downgoing nuclearites 
with /3 = 2 X 10"'^ at ground level. Nuclearites of such velocity could have 
galactic or extragalactic origin. The MACRO direct limit is shown along 
with the limits from Refs. [13] ( "Nakamura" ) , [14] ("Orito") and the indirect 
mica limits of Refs. [15,16] (dashed line). The MACRO limit for nuclearite 
masses larger than 5 x lO^^GeV/c^ has been extended and corresponds to 
an isotropic flux. 
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